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The low temperature phase diagram of ID weakly disordered quantum systems like charge or 
spin density waves and Luttinger liquids is studied by a full finite temperature renormalization 
group (RG) calculation. For vanishing quantum fluctuations this approach is amended by an exact 
solution in the case of strong disorder and by a mapping onto the Burgers eguation with noise 
in the case of weak disorder, respectively. At zero temperature we reproduce the quantum phase 
transition between a pinned (localized) and an unpinned (delocalized) phase for weak and strong 
quantum fluctuations, respectively, as found previously by Fukuyama or Giamarchi and Schulz. 

At finite temperatures the localization transition is suppressed: the random potential is wiped 
out by thermal fluctuations on length scales larger than the thermal de Broglie wave length of the 
phason excitations. The existence of a zero temperature transition is reflected in a rich cross-over 
phase diagram of the correlation functions. In particular we find four different scaling regions: a 
classical disordered, a quantum disordered, a quantum critical and a thermal region. The results 
can be transferred directly to the discussion of the influence of disorder in superfluids. Finally we 
extend the RG calculation to the treatment of a commensurate lattice potential. Applications to 
related systems are discussed as well. 
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I. INTRODUCTION 

The collective behavior of condensed modu- 
lated structures like charge or spin density waves 
(CDWs/SDWs)i2£, flux line lattices^ and Wigner 
crystals^ in random environments has been the subject of 
detailed investigations since the early 1970s. These were 
motivated by the drastic influence of disorder: without 
pinning CDWs would be ideal superconductors whereas 
type-II superconductors would show finite resistivity. In 
three dimensional systems the low temperature phase 
of these structures is determined by a zero temperature 
disorder fixed point resulting in quasi-long-range order 
and glassy dynamics (for recent reviews and further 
references see 00). In two dimensions this fixed point 
is extended to a fixed line which terminates at the glass 
transition temperature^. In the low temperature phase, 
correlations of the positional order decay slightly faster 
than a power law and the linear resistivity vanishes^. 

In one dimension the situation is different: the glass 
temperature is shifted to T — 0. Nevertheless, there re- 
mains a residual trace of disorder which is reflected in the 
low temperature behavior of spatial correlations and the 
dynamics&S. Clearly, at low temperatures also quantum 
fluctuations have to be taken into account. Disorder and 
quantum fluctuations in ID CDWs at zero temper ature 
have been considered previously (see, e.g., |lOlll| ) and 
an unpinning (derealization) transition as a function of 
the strength of quantum fluctuations was found. Finite 
temperature effects were partially incorporated by trun- 
cating the renormalization group (RG) flow at the de 
Broglie wave length of the phason excitationsii. How- 
ever, for a complete study of the thermal to quantum 



crossover, quantum and thermal fluctuations have to be 
considered on an equal footing^, which is the main aim 
of this paper. 

Experimentally, quasi-lD behavior can be seen in real 
materials, e.g., in whiskersi^, hairlikc single crystal fibers 
like NbSe3, with a transverse extension smaller than the 
correlation length or in chain like crystals with weak in- 
terchain coupling. In the latter case there is a large 
crossover length scale up to which ID behavior can be 
observed^. 

The results obtained for the CDWs or SDWs have a 
large number of further applications on disordered quan- 
tum systems: they relate, e.g., to the localization transi- 
tion of Luttinger liquidsiS*ii, superfluids^, tunnel junc- 
tion chains^, Josephson coupled chains of these systems, 
if the coupling is treated in mean-field theorjiifi, and 
CDWs in a lattice potential. However, we will use the 
terminology of CDWs in most parts of this paper. 

The remaining part of the paper is organized as follows: 
In section[n]we give a detailed introduction to our model 
and the notation used in this paper. We also briefly dis- 
cuss the influence of Coulomb interaction on the proper- 
ties of the system. In section lTTTl the influence of the disor- 
der is studied in detail. Using an anisotropic momentum 
shell renormalization group calculation, in which the full 
Matsubara sum over frcquncies is performed, we obtain 
flow equations for the effective strength of the disorder, 
thermal and quantum fluctuations (i.e., the interaction 
strength in the case of Luttinger liquids). These are dis- 
cussed first in the case of zero temperature and agree- 
ment with previous results is obtainecU2*ii. At finite 
temperatures the disorder always renormalizes to zero. 
In the classical limit two more methods are applied: (i) 
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at low temperatures and strong disorder the ground state 
of the model is calculated exactly, (ii) For weak disorder 
and strong thermal fluctuations a second RG calculation, 
based on the mapping onto the Burgers equation with 
noise, is applied. Using all these findings, the phase dia- 
gram of the density-density correlation function is stud- 
ied in section llVl The main result of this section is the 
calculation of the low temperature quantum crossover di- 
agram for one-dimensional CDWs. In the following sec- 
tion [V] we discuss briefly the application of the results 
to superfluids by using the mapping to CDWs. Some of 
these results were previously presented in Ref. Il6l 

The influence of a commensurate lattice potential on 
a free density wave is considered in section IVll The full 
finite temperature renormalization group flow equation 
for this sine-Gordon type model are derived and resulting 
phase diagram is discussed. 

In the appendices we present the calculation of the 
renormalization group flow equations and the derivation 
of the correlation function in the strong and weak pinning 
limit in some detail for the interested reader. In the final 
section of the appendix we list all symbols used in this 
paper with corresponding references in the paper. 



II. MODEL 
A. charge and spin density 

In this section we derive the effective Hamiltonian 
which will be the starting point for our further treatment. 
The strategy of the calculation is therefore separated into 
two steps. In the first step the system is treated in a 
mean-field (MF) type approximation applied to a micro- 
scopic Hamiltonian. This leaves us with a slowly vary- 
ing complex order parameter field for which we derive 
an effective Hamiltonian. The second step involves the 
consideration of the fluctuations of this order parameter, 
which is the topic of this paper. 

We briefly summarize now the result of the mean-field 
calculation: Well below the mean-field condensation tem- 
perature Tmf of the CDW, the underlying lattice will be 
periodically distorted with a period A which is related to 
the Fermi wave vector kp by A = n/kp. This distortion 
of the lattice leads to the formation of a gap in the dis- 
persion relation at k = ±kp which is (in one dimension) 
proportional to the amplitude of the lattice modulation. 
For small displacements (which are typically smaller than 
1% of the interatomic spacing^), the increase of the elas- 
tic energy is smaller than the gain of electronic energy 
due to the formation of the gap and hence an instability 
is favored. The period of the CDW depends on the band 
filling factor (via kp = n / X) and is in general at arbitrary 
band filling incommensurate with the undistorted lattice 
(with lattice constant a). 

In (quasi-) one-dimensional systems^ also SDWs can 
be found, but in contrast to CDWs they arise due to 
electron-electron and not to electron-phonon interaction. 



A SDW can be considered to consist of two CDWs, one 
for spin-up and another for spin-down electrons (see, 
e.g., Fig. 5 in Q). Therefore the spatial modulation of 
SDWs is characterized by a wave vector Q = 2k f, as for 
CDWs. 

The charge or spin density p(x, t) can be written in the 



p(x,t) = (1 + Q 1 d x ip(x 1 t))[p + pi cos(ptp(x,t) +pQx)} 

(1) 

where po = Qf(T)/ir and p\ = 2| A\/(irgvp ). g is the 
dimcnsionless electron-phonon coupling constant and vp 
the Fermi velocity. p\ is proportional to |A|, the CDW 
gap or the amplitude of the complex (mean field) order 
parameter 



A(x,t) = \A{x,t)\e^ x ' l \ 



(2) 



/ is the condensate density (Y = 1 
Yoshida function) related to |A| by 



/ is the so-called 



f(T) 



ttT 



E 



(|A|A) 2 



\2vV2' 



2imT/h. 



(3) 

The condensate density approaches 1 for T — > and 
f(T) ~ 2(1 - T/T C MF ) ~ |A| 2 for T -> T C MF . T C MF 
denotes the mean-field transition temperature. For quasi 
one-dimensional systems p\ has an additional factor £~ 2 
(the area perpendicular to the chain) . 

Note, that is correct for arbitrary band filling and, 
to be more precisely, is the particle density of the charge 
or spin carrying particles. The particle current density j 
follows from J3J) as j = —pgip/Q. 

Because Akp modulations of SDWs or CDWs are also 
possibloi^, we introduce the factor p in the argument 
of the modulating cosine function, i.e., for CDWs and 
SDWs p is usually 1, but can also be 2 or greater. 



B. Hamiltonian 

In the following we use a minimal model for the low 
energy, long wave length excitations of the condensed 
charge density wave. Since fluctuations in the amplitude 
|A| are suppressed, because these are massive, we take 
into account only fluctuations in the phase ip. Clearly, 
such an approach breaks down sufficiently close to the 
mean-field transition temperature T c Mi? . Neglecting fluc- 
tuations in |A|, the Hamiltonian of the CDW is given by 



(4) 
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with 



T~to = 



7~Lu 



H w = 



dx\ 



- P 2 + (d x 0) 2 , (4a) 
cJ 

dxU(x)p(x) , U(x) = UiS(x - Xi){4b) 

L 

— J dxW cos (qtp(x)) . (4c) 





Ho describes the phason excitations of the CDW, where 



4^/(T) denotes the elastic constant. 



vf/ 'y/T+ (2|A|7 hhjpq) 2 J (gf) is the effective velocity of 
the phason excitations with oj p q the phonon frequency. 
For CDWs (2\A\/hu pQ ) 2 /(gf) > 1 is typically fulfilled 
and hence quantum fluctuations are weak. 

P is the momentum operator, corresponding to 
the phase <p, with the standard commutation relation 
[P(x),0(x')) = i6(x-x') 

Tin results from the effects of impurities with random 
potential strength Ui and positions xi. The potential 
strength is characterized by U< 



and UiUi 



U 2 6- ■ 



and includes a forward and a backward scattering term 
proportional to po and pi, respectively. The disorder 
average of the impurity potential U(x) follows then to be 
given by U{x) = and 



U(x)U(y) 



-6(x - y) 



(5) 



imp 



We will further assume, that the mean impurity dis- 
tance limp = L/Ni mp is large compared with the wave 
length of the CDW and, in most parts of the paper, that 
the disorder is weak, i.e., 



1 < kmpQ < cQ/{U lmp pi) 
In this case the Fukuyama-Lee length 

2/3 



L 



FL 



\U imp pip 2 



(6) 



(7) 



is large compared to the impurity distance li mp . 

The third term in (@J, Ti W: includes the influence of a 
harmonic lattice potential. This term will be discussed 
section IVll in greater detail. 

Our model Q includes the four dimcnsionless param- 
eters 



t 
K 



T/ttAc , 
hv/wc, 

{UimpPl) 

W/ttcK 2 , 



(8a) 
(8b) 

(8c) 

(8d) 



which measure the strength of the thermal (t) , quantum 
(K) and disorder fluctuations (u) and the periodic po- 
tential (w), respectively. A = ir/a is a momentum cut- 
off. The classical region of the model is given by K <C t 
which can be rewritten as the condition, that the thermal 
dc Broglie wave length 

A T = hf3v = K/(tA) (9) 

of the phason excitations is small compared to a. 

At T = 0, AT-values of the order fO" 2 to ICT 1 and 
1, have been discussed for CDWs and SDWs, respec- 
tivelyiSS. It has to be noted however, that the the ex- 
pressions relating c and v to the microscopic (mean-ficld- 
like) theory lead to the conclusion that K and t diverge 
by approaching T C MF , whereas the ratio K/t remains fi- 
nite. 



C. Coulomb Interaction 

We could also add a Coulomb interaction term to our 
model Q which can be written as 



^4 

2 



dx I dx' p(x)V c (x — x')p(x') . 



(10) 



where V c is the Coulomb potential. In all dimensions 
the unscreened potential has the form e 2 /r. If we as- 
sume, that the quasi one-dimensional system has the fi- 
nite width £, V c can be written as^iSS 



V c °(x) 



1 



v/x 2 + C 2 



^e^V^fc) with (II) 



L ^ 

k 



K (fc) = 2e 2 K (ICfc|), 



(12) 



where Kq is a Bessel function with Kq(x) ~ — ln(x) for 
x < 1. 

In general the Coulomb potential is screened and can 
be written asS 



V c {k,u) = 



V c °(k) 



i + v?(*)n(fc,w), 



(13) 



with the momentum and frequency dependent polariza- 
tion operator U(k, u) = (p(0, 0)p(k, u)) . 

If we only consider the static case ui = we can distin- 
guish two limiting cases: First, if the typical range A eJ y of 
the screened Coulomb potential V c is much smaller than 
the mean electron distance, the potential can be assumed 
to be a delta distribution and H c can be approximated 
by 



% c w ~y I dx (—d x tp(x) 



(14) 



with x = \ J dxV c (x). The cos-terms (. . .) from the den- 
sity can be neglected due to strong fluctuations. There- 
fore the Coulomb interaction gives only an additional 
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contribution to the elastic constant of the initial model: 
c = A^e./ + |j. For x > the Coulomb interaction is 
repulsive, which leads to an increase of c and therefore 
a decrease of the dimensionless parameter K, i.e., the 
quantum fluctuations will be reduced by the Coulomb 
interaction. In the case \ < (attraction), K will be 
increased. Keeping this consideration in mind, we will 
not further include TL C in the model explicitly. 

In the other case - with weak screening - V c (k) m 
V c °(fc) shows the dispersion given in (|12fl and in gen- 
eral, the details of the fc-dependence depend not only on 
the transverse extension £ of the quasi one-dimensional 
system under consideration but also on the screening 
length^,. 

However, the logarithmic fc-dependence will only 
weakly affect our RG-analysis, but may suppress phase 
transitions, as discussed later in section ITlI Bl 

Coulomb interaction is also important if one consid- 
ers multi-channel systems^ or the effect of the non- 
condensated normal electrons. 



which starts from a path integral formulation of the par- 
tition function corresponding to the Hamiltonian We 
begin with the renormalization of the disorder term and 
put w = in the following. The system is transformed 
into a translational invariant problem using the replica 
method, in which the disorder averaged free energy is 
calculated, using 



T = -TlnTre-' s / s = — T lim - fire- 5 '"' V fi - l) , 



n— »o n 



(15) 

which defines the replicated action S^ n \ is given by 



III. RENORMALIZATION GROUP 
TREATMENT OF DISORDER 

A. Flow equations 

In order to determine the phase diagram we adopt a 
standard Wilson-type renormalization group calculation, 



{t))H u [pu{t')] 



(16) 

where £q is the Lagrangian corresponding to Ho, J = 

J Q h ^ dr and a, v are replica indices. Using (0) and con- 
sequently neglecting higher harmonics (2pQ-modes) one 
finds 



Hu[<Pa(T)YH u [<Pv(T')] = 2™ PPl j dx jcOSf^a^, r ) ~ fv(x,r')J + ^° 2 d x (p a (x,T)d x ip v (x,T')\ . (17) 

imp ^ 

Together with (|16f) one obtains the following form 
h ~ 2ttK 



LA K/t . K/t 

J dx J dT \ {dxfa) 2 + {d r f a ) 2 S av - J dr' u 2 cos p (p a (x, t) - <p u (x, t' 



+ad x ip a (x, r)d x tp v (x, t') 



(18) 



with a = 2u 2 (p A/piQf 



r 



Note, that we introduced dimensionless spatial and imag- 
inary time variables, 

Ax — > x , 

Al!T — > T , 

which will be used throughout the paper - beginning 
here. Furthermore all lengths (e.g. correlation lengths, 
At, Lfl, li mp and L), wave vectors (e.g. k, kp and 



Q) and Matsubara frequencies are dimensionless accord- 
ingly, from now on. Additionally we rescale the elastic 
constant 

Ac — > c , 

for convenience to avoid the appearance of A. 

Integrating over the high momentum modes of ip(x, r) 
in a momentum shell of infinitesimal width l/b < \q\ < 1 
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but arbitrary frequencies and rescaling x 



x/b, 



r/b, we obtain the following renormalization 



group flow equations (up to one loop): 
dt 



dl 
dK 

~dT 
df_ 

~~dT 

(1(7 



l pVKB J P 'K,§ 



3 — — — coth — 
2 2t 



coth 



K 
2? 



(19a) 
(19b) 

(19c) 
(19d) 



where I = In b. For details on the RG calculation we 
relegate to Appendix lAl where we have written the RG- 
flow also for dimensions < \d — 1| -C 1. 

For legibility we have introduced the following func- 
tions: 



Bi{v,y) = 



dr 
o o 



, x 9i(T,x) cosh (y-r)^ (2Q) 



T(t,x) 



with 



T(t, x) coshy 

( U\ z ( ' KX 7TT 

1 + I — cosh cos — 

Vtt/ \ y y 

9o(t,x) = S(x)t 2 , 



-i i//4 



(21) 



K- 
It > 



for K -> (see Fig. in 



Note, that B (p 2 K, ■ 
Appendix [X| . 

There is no renormalization of t (i.e., of the elastic con- 
stant c) because of a statistical tilt symmetry^. Note 
that l|18f) is written in rescaled dimensionless parameters 
and the different renormalization of the kinetic and elas- 
tic term is reflected in the different renormalization of v 
and c, i.e., K and t, respectively. 

From the flow equation for u 2 l|19cfl one sees directly 
that, depending on the sign of the prefactor, the behavior 
changes from increase for small t and K to decrease for 
high K or t. 

There is no first order RG correction to a and the 
change of a with length scale is simply given by rescal- 
ing, see (|19d|l . The two-loop contribution to a is much 
more involved than the one-loop contributions for the 
other flow equations and gives no qualitatively different 
result for the flow of a. As seen from (I19dll . the forward 
scattering amplitude always increases as cr^e 1 on larger 
length scales and is therefore not well controlled in the 
RG sense. But, since the flow of a docs not feed back 
into the other flow equations it has only minor relevance 
for our considerations. And indeed, we can get rid of the 
forward scattering term / /ttII(x)^ by introducing the 
field (pb{x) by 



<p{x) = <p b (x) - ip f (x) , Lp f {x) 



dyc(y) , (22) 



where c(x) = }cA > w ith c(x) = and c(x)c(x') = 
ja5(x — x'). This can be easily verified by inserting 



this decomposition of (p(x) into the initial Hamiltonian 
PJ) written in dimensionless units, and using (jSJ and the 
definition of a for deriving the averages of c(x). Note, 
that x is dimensionless. 



TABLE I: Notation guide, 
symbol here G. fc S.", Ref. [Hi Haldane, Ref. [311 



P 
K 

v 
c 

V 



V2<j>, V2/nd T (f> 



2K a 



2tvK p 
1 



■ npox, TT 



tlVN /V 

2 



"charge operators 



B. Zero temperature - a review 

The special case t = was previously considered, e.g., 
m E3 (for a better comparison see our notation guide 
given in table QJ. 

The flow equations for K and u at zero temperature 
read: 



with 



dK 
~dT 

dl 



Bo(v, oo) 



-A 2 a^ (p 2 a;oo). 

o P 2 K 



dTT 2 e- T 



[l + r 2 /2] 



-u/i 



(23a) 
(23b) 

(24) 



The corresponding flow equation for K obtained in [Tlj 
deviates slightly from (|23a|l , which can be traced back to 
the different RG procedures. In ^lj the authors per- 
formed the RG at strictly zero temperature and used 
a symmetric, circular shape of the "momentum-shell", 
i.e., treated the model as effectively isotropic in the 1+1- 
dimcnsional space-time. 

This procedure may be a good approximation at zero 
temperature, but if one considerers finite temperatures 
this does not hold anymore, since the extension in in- 
direction is now finite. As a result, there is a region 
tt I L < \k\ < tt/\t where fluctuations are mainly one- 
dimensional and purely thermal. This region was disre- 
garded in previous treatments. As we will see, fluctua- 
tions from this region have an important influence on the 
overall phase diagram. 

The critical behavior is, however, the same: there is a 
Kosterlitz-Thoulcss (KT) transition at the phase bound- 
ary K u between a disorder dominated, pinned and a 
free, unpinned phase which terminates in the fixed point 
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K* = 6/p 2 . One can derive an implicit equation for K u 
by combining (|23a[) and (|23b|) to a differential equation 



du 



dK p 2 r]K 
which has the solution 



(K - K) . 



_K* u f K-K 
P 2 V V K l K 



\K)-ul = 



(25) 



(26) 



where uq and Kq denote the bare values of the dis- 
order and quantum fluctuation, respectively, and -q = 
Bq(p 2 K*,oo). Then, K u is implicitly given by 

u 2 (K u ) = ^(^^-^) , (27) 
P 2 V \ K* K*J 

where the initial condition u 2 (Kq = K*) = Uq = is 
used. The KT-fiow equations at K* can be recovered by 
defining 

P 2 K „ 



2 X 2 = 2 p4??u2 



with |7| < 1 . This yields 

d'y 

U 

d X 2 



(28a) 
(28b) 



which are exactly the flow equations obtained by Koster- 
litz and Thouless^. 




FIG. 1: Schematic zero temperature phase diagram in d = 1 
and close to d = 1 dimensions (see text), u and K denote the 
strength of the disorder and quantum fluctuations, respec- 
tively. 



Fig- E>> whereas for e > 0, K and u always flow to 
the strong pinning fixed point (at K = and u — > 00; 
right inset), i.e., quantum fluctuations are too weak to 
renormalize the random potential to zero. The zero tem- 
perature transition disappears therefore for d > 1, since 
the fixed point lies in the unphysical u < region of the 
K-u parameter space. This can easily be verified by us- 
ing the rescaling of K given in eq. (|AI2c|l of Appendix 1X1 
which results in the flow equation (|A16(1 . In general this 
discussion applies to the localization transition as well as 
to the Mott transition (see discussion of the the lattice 
potential). Note, that the flow for d ^ I is qualitatively 
different from that discussed in Ref. [23 ■ 

If one includes the effect of Coulomb interaction in 
d = 1 dimension, phase fluctuations of the free phase 
field increase only as (T=0) 



(Mi,0)-#0)) 2 )~ Kln^ 2 \x\ 



(29) 



As a result, phase fluctuations are too weak to suppress 
the disorder even for large values of K and the system 
is always in the pinned phase. The phase diagram is 
therefore similar to that in d > 1 dimensions. 

In the pinned phase the parameters K and u flow into 
the classical, strong disorder region: K — > 0, u — > 00. 

Integration of the flow equations gives for small initial 
disorder and K K* an effective correlation length or 
localization length 

(l-K/K'J- 1 



FL 



(30) 



at which u becomes of the order unity. This can be ex- 
tracted from (|23b(l neglecting the flow of K. 

A better approximation of £ u which takes also the flow 
of K into account can be obtained by replacing u 2 in 
the flow equation for K Ij23a(l by the expression given in 
l|2B). We still use the approximation, that K deviates 
not much from the bare value Kq which is the 
long as UqI f . Then, the solution for K{1) is given by 



( 1 - y«o^ 



K{1) 

which yields a solution of l|23bl 
u 2 (l) 



(31) 



In 



3 - ^K a ) I + ^r,K u 2 l 2 . (32) 



With u 2 (\n(£ u )) sa 1 , the correlation length £ u is defined 
by 

= ln«g+ U - ^K \ ln(e u ) + y^o«o(hi(£u)) 2 (33) 



Under the assumption, that a small deviation from the 
dimension d = 1 changes only the naive scaling dimen- 
sions of the fields, our results can be extended also to 
d = 1 + e dimensions (For details see Appendix B). The 
zero temperature phase diagram is modified and illus- 
trated in Fig. [I] For e < the fixed point at (K = K*, 
u = 0) is shifted to positive u-values (see left inset of 



which yields 



\J a 2 — 4b In it§ — a 



2b 

In Uq 

_ p 2 K 

2 



,2\2 



V 9 (In tin) 

^Kml - K ° _ - 3 ,(34) 
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where the first term of the r.h.s. gives the result (|30f) . 

Close to the transition line £ u shows KT behavior. For 
K > K u , £ u diverges and C(x,t) ~ K(l — ln|z|)ln|z| 



where \z\ = \Jx 2 + r 2 (cf. section ITVl 



C. Strong pinning limit: Exact ground state 

For large values of u our flow equations break down. 
Qualitatively the flow is towards large u and small K. 
We can, however, find the asymptotic behavior in this 
phase by solving the initial model in the strong pinning 
limit exactly. To find this solution we will assume strong 
pinning centers and weak thermal fluctuations: 



which defines the exact ground state of the classical 
model: If we use l|38(l one finds for the optimal value 
of the ni 's 



n Q i+1 = n? 



pQ 

2tt 



{c-i -\- limp] 



which leads, using (|37J) . to the exact classical ground 
state 



pQ 

2tt 



\£j ~t~ limp) 



where n\ has an arbitrary integer value. 



(43) 



oo and c/Ump S> T. 



(35) 



To treat this case we go back to the initial Hamilto- 
nian (0J (with W = and the kinetic term also vanishes 
because of K — > 0). For strong disorder it is convenient 
to integrate out the phase field f{x) at all points which 
are not affected by the impurities. Then the effective 
Hamiltonian takes the fornA 



n 



N 

E 



2 



+ Uip(xi) > , if 



Under condition (|35|l . ipi only takes values obeying 



(36) 



p{<Pi + Q x i) = 2^77,^ + 7r with rij S Z integer (37) 

which minimizes the backward scattering term. Defining 
hi and q by 



n i+ i 



= hi 



pQI 



imp 



2n 



l imp ~ 



(38) 



with < x\ < X2 < ■■■ < xn+i < L, the effective 
Hamiltonian can be rewritten as 



(39) 



Here [x] denotes the closest integer to x ( Gaussian brack- 
ets) : 



[x] = ni for x € 



1 1 

m , m H — 

2' 2 



, to € Z 



and 



imp 



27T 



2- 



(40) 



(41) 



such that |7| < 5. 

Because thermal fluctuations are small compared to 
the elastic energy, see <|55|) . (hi — ^^li _ -y) takes on its 
minimal value, which is given by 

pQti , 



2tt 



(42) 



D. Finite temperature and crossover diagram 




FIG. 2: Typical flow diagram for w = in the three dimen- 
sional parameter space of K, u and t, proportional to the 
strength of quantum, disorder and thermal fluctuations, re- 
spectively. 

At finite temperatures thermal fluctuations wipe out 
the random potential which lead to the pinning of the 
CDW at t = and K < K u . Thus there is no phase 
transition any more, in agreement with the Landau the- 
orem. The system is always in its delocalized phase even 
if the disorder may still play a significant role on inter- 
mediate length scales. 

In the special case K — > the flow equation (|19c|) 



reduces to 



u\l) 



p 2 t u 2 with solution 



u 2 e 3I-p 2 t (e , -l) 



(44) 



If we write t — toe 1 we may express I by t and hence we 
may write u 2 as t-dependent function: 



u\t.) = ul{t/t fe-v 2 ^ 



(45) 



which is plotted in Fig. [21 in the t-u plane. 

One sees that the flow of the disorder has a maxi- 
mum at t = 3/p 2 or / = \n(3/(p 2 t a )), if t a < 3/p 2 . 
For finite K, the RG flow of u in the region K < K u 
first increases and then decreases. The region of in- 
crease in the K-t plane is implicitly defined by M u = 
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{(K,t)\K* > ifcothf > 0}, i.e., the positions of the 
maxima of u 2 [K,t) are located on the boundary of M. u 
defined by K* — K coth || . 

The correlation length £ can be found approximately 
by integrating the flow equations until the maximum of 
u(l) &ndt(l)/(l + K(l)) is of the order one (see discussion 
in section flvl . This can be done in full generality only 
numerically (see Fig. |3J). 



t A 



'u 
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0.02 
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M 

1 classical 
, disordered 
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quantum 
critical 



quantum 
disordered 
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1 1-2 1.4 K" 



FIG. 3: The low temperature crossover diagram of a one- 
dimensional CDW. t and K are proportional to the temper- 
ature and the strength of quantum fluctuations, respectively. 
The amount of disorder corresponds to a reduced tempera- 
ture t u ~ 0.1. In the classical and quantum disordered re- 
gion, respectively, essentially the t = behavior is seen. The 
straight dashed line separating them corresponds to At ~ 1, 
i.e., K w t. In the quantum critical region the correlation 
length is given by At- Pinning (localization) occurs only for 
t = 0,K < KZ. 

It is however possible to discuss several special cases 
analytically. The zero temperature correlation length can 
still be observed as long as this is smaller than the ther- 
mal de Brogle wave length At which can be rewritten 
for K not too close to K* as t < tx ~ Ktu K I K ^ 
with t u « L~j? L , where we defined tx via £ u = anal- 
ogously to the definition of At, and used (|30|l . We call 
this domain the quantum disordered region. 

For K > K u the correlation length £ is given by At 
which is larger than given by purely thermal fluctuations. 
For scales smaller than At, the phase correlation function 
still increases as ~ In \z\ with a continuously varying co- 
efficient K e ff(uo), as will be discussed in detail in the 
next section. In this sense one observes quantum criti- 
cal behavior in that region, despite of the fact, that the 
correlation length is now finite for all values of K [T3| . 

In the classical disordered region tx < t < t u the cor- 
relation length is roughly given by Lfl as follows from 
previous studies— or by solving u 2 (ln(£)) ~ 1 using 



-2/3 



-v to 



-2/3 



(|44|1 for small to yielding £ s= 
L FL (np 4 ) 1/3 . Note, that t K w K for small K. 

In the remaining region t>t u , the thermal region, we 
apply the mapping onto the Burgers equation (see sec- 
tion ■ I n this case the RG-procedure applied to this 
equation becomes trivial since there is only a contribu- 
tion from a single momentum shell and one finds for the 
correlation length w § /(T)t[l + l/2[t tl /( 7 rp 2 i)] 3 ]A. 

The phase diagram depicted in Fig. [3] is the result of 



the numerical integration of our flow equations and shows 
indeed the various crossovers discussed before. 

In the high temperature region (t ^> K) the flow equa- 
tions can be solved explicitly. For u 2 (l) we get the same 
result as given in (|44|l and the flow equation for K re- 
duces to 



dK 

~dT 



2 U {2tf 



(46) 



where we used B (p 2 K, § -> 0) = (K/2t) 4 . The solution 
of this equation is given by 

3p*u§ 



K(l) 



ml 



e p to Ei(p 2 t ,p 2 t e l 



-1/3 



(47) 



with the incomplete exponential integral function Ei(a, b) 
defined by 



Ei(a,6) = / dte~ f 

J a 



/*■ 



One observes that K(V) saturates very quickly to the 
value K(oo) < K - 



IV. CORRELATION FUNCTIONS 

In this section we discuss the density-density and the 
phase correlation functions in more detail and summarize 
all correlation lengths in the various regimes - partly al- 
ready used in the last two sections. 

The (full) density-density correlation function is de- 
fined as 



S(x,t) = (p(x,r)p(0,0)) 



(48) 



where p{x,r) is given in In the following we re- 
strict our considerations to the ( charge ) density wave or- 
der part of 5, which is the term proportional to p\, i.e. 



S 1 (x, t) = p\ (coap(ip(x, t) + Qx) cospip(0, 0)) , (49) 

which defines the type of order of the density wave: If 
it decays algebraically we have quasi long-range order 
(QLRO), an exponential decay over a correlation length 
£ corresponds to short-range order (SRO). The omitted 
parts of S decay faster than Si j2^ ■ 
5i can be rewritten as 

2 

S^x, T ) = ^ ( e lpQx ^p(v(*,t)- v (o,o))^ + 

e -ipQx / e -ip(<p(x,T)-<p(0,0))^ j ; (5Q) 

and under the assumption that we can use a gaussian 
approximation for the averages, we obtain 

Si(x, t) ~ p\ coa{pQx)e-^( {lp(x ' T) - lp{0 ' 0))2 ) . (51) 
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From now on wc focus on the phase correlation func- 
tion 



C(£,t) = ((^,t)-^(0,0)) 2 ) 



Note, that with this result we have neglected the al- 
gebraic decay for small x < At- Therefore a bet- 
ter interpolation formula for the correlation length is 
(52) £ sa -4(£r + At), which takes the slow decay for small x 



and discuss it in various limits. Combining (|51|l and Ij52(l 

we can extract a correlation length from the relation 



r 1 = lim %-C(x,0) 



Disorder-free case 



(53) 



We start with the most simple case u = 0. Then, 
the correlation function in dimensionless units follows di- 
rectly from the action Sq written in momentum space: 



C (X,T) = — ^ 2 
fc,n " 



(54) 



with Matsubara frequencies aj„ = 2 7m/ At and momenta 
& = &m = 2Ttm/L). 

The sums over n and k (i.e., m) can be performed 
approximately for sufficiently large x and r and one ob- 
tains^ 



C (.t,t) 




(55) 



/2ttt 
I At" 



The behavior of this function is considered in the fol- 
lowing cases: 

(i) At zero temperature (At — > oo) (|55|l reduces to 



C ^r)^^ln(i[x 2 +T 2 ]+1 



(56) 



i.e., the correlation function has a logarithmic depen- 
dency on x and r and leads to an algebraic decay of 
Si, i.e., the system shows QLRO. 

(ii) At finite temperatures we can distinguish between 
length scales smaller and larger than At- 

In the first case x <C At and r <C At the cosh and cos 
term can be expanded to second order in the arguments 
and one gets the same logarithmic function as in the zero 
temperature case. In the opposite case x At , which is 
the usual case at high temperatures, the cosh term can be 
approximated by the exponential function and one finds 
a linear dependency on x: 



Cq{x) w ntx = Tx/c 



p 2 irt 
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(57) 



i.e., 5*i decays exponentially (SRO) over a characteristic 
length £ ~ t~ l . The same result is obtained for the limit 
K — > at a fixed, finite temperature. 



into account. In terms of the length-scale dependent t(l) 
this rewrites to 



t(l = ln(£)) — K + 1 



(58) 



i.e., the correlation length is reached, if t(l)/(l + K) is of 
order one. 

The change from QLRO on small length scales x < £ 
to SRO on large length scales becomes clear if one consid- 
ers the cylindric topology of the system in space-time at 
finite temperatures: As soon as one reaches length scales 
of the order of the perimeter of the cylinder, which is 
At, starting from small scales, the system changes from 
two-dimensional to effectively one-dimensional behavior. 



B. 



Finite disorder 



If u is finite the action of the system has a forward 
and a backward scattering part. With the decomposi- 
tion (|22|1 . the phase correlation function divides into two 
parts: 



C(x,T) = C b (x,T)+C f (x) 



(59) 



and has therefore always a contribution Cf(x) ~ |s|/£/ 
with £7 ~ o~(l = ln|a;|), i.e., the density wave order 
has always an exponentially decaying contribution and 
we can define 



Si(x,r) = f p (x)e~ 



-C b (x,T) 



(60) 



with f p (x) = p\ cos(pQx)e~ E ^^^ f . However, since 
Cf(x) is not r-dependent, it will not influence the dy- 
namical properties of the system. Therefore all further 
remarks about phase correlations refer to Cb{x,r) and 
consequently we will drop the subscript b in the follow- 
ing. Again we examine the T = and finite temperature 
cases: 

(i) At zero temperature we have to distinguish between 
three i^T-regimcs: For K > K u the disorder becomes ir- 
relevant under the RG flow and we can use the zero tem- 
perature, disorder free result for the correlation function 
with the pre-factor K replaced by an effective quantity 
K e ff(l = In 2) on a length scale z = \J x 2 + t' 2 , defined 
by the flow equation for K. This effective K saturates 
on large scales at a value K e ff(uo), which may be seen in 
Fig. |3 Therefore we have QLRO in this K region. 

For < K < K u wc integrate the flow of u until it 
reaches a value of order one, starting at small uq, which 
defines the localization length £ M (see section UlI B|) . i.e., 
the correlation function behaves like C(x,t) ~ 
i.e., we have an additional (to Cf) exponentially decaying 
contribution to Si . 
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C. Strong disorder 



and for l imp » Q 1 , i.e., a«l: 



In the last region K = we come back to the strong 
pinning case, discussed in section IIII CI before and cal- 
culate the pair correlation function exactly. Taking into 
account that the hi's arc independent on different lat- 
tice sites, i.e., hihj <x 8ij, the (discrete) phase correlation 
function is given by 



71+ 1 




where we used (|42|) for the second equality. For eval- 
uating the disorder average in this expression, one has 
to take into account the order statistics of the impurity 
distances e^. In the thermodynamic limit the probability 
density function for the e,'s can be rewritten as 



-l 

imp — 1 7 



/ \ imn 

p{ei) « e 

e 



-kmp < e, < oo . (61) 



Then, the correlation function can be explicitly written 

as 



2 \ 4tH 



dxe x 

\2a 



where we introduced the parameter a 



x 
2a 



(62) 
and sub- 



pQlimp 

stituted x = l~^ p £i + 1- This integral can be evaluated 
exactly, which leads to the following exact expression for 
the pair correlation function at zero temperature, written 
in a continuum version: 



C(x,t) = — ( 
pa 



sinha 



p 2 £ 



pQ^imp 

(63) 

A more detailed derivation of this result is given in Ap- 
pendix \n\ 

Finally we want give an interpolating expression for 
C(x, t) from T = to high temperatures T ^> cj (hmpP 2 ) 
starting with the result (|63|l . In the latter case we may 
neglect the discreteness of hi and hence 



(tp n+1 - (pi) 
4tt 2 " 



p 2 i 

imp 

T 

— \x\ = irt\x\ , 
c 



— — — — In / dhe-Zi x - h2 



(64) 



(mx) - ^(o)) 2 ) « (£i-± + 



3p 2 imp 18Qp 4 Q2 

and hence the correlation length acquires the form 



1 - 



sinh(a) 



(66) 



(67) 



Note that U m pQ > 1, he., a < ir/p and £t 3> kmp- An 
approximate crossover to the weak pinning limit follows 
by choosing k mp ~L FL . 

(ii) At finite temperatures the parameter K satu- 
rates at an effective value K e -g{uo) on large length scale. 
Therefore the correlation function for small disorder is 
given by (|55|l with K replaced by K{1 = lnz). 

In the region M. u of the K-t plane (see section Ull PL 
u still increases and we can find an effective correlation 
length by comparing the length scales on which u(l) or 
t(l)/(l + K(l) become of order one. Then, the correlation 
length is the smaller length of these two. 



D. Burgers equation 

For K = 0, high temperatures but weak disorder we 
adopt an alternative method by mapping the (classical) 
one-dimensional problem onto the Burgers equation with 
noised. With this approach one can derive an effective 
correlation length given by 



1, £r , 
2 [ 2L FL ' 



(68) 



where £r <C Lfl, which changes the prefactor of the free 
correlation function at high temperatures 1)57(1. The full 
calculation of this result can be found in Appendix IH1 



TABLE II: Overview of the dimensionless correlation lengths. 



length 


description 


eq. 




weak pinning/high temp, length 




£/ 


forward scattering length 






strong pinning length 






high temp. /disorder free length 






disorder localization length 






lattice pot. correlation length 


sec. EH 



With Xi = Tptfkmp+ei) ■ 

A plausible interpolation formula is then given by 

Q 



Mx) - <p(o)y 



2Q 2 l 



imp 



l 



sinh(c^ 



T 
— I 



(65) 



V. SUPERFLUIDS 

Next we consider the application of the results ob- 
tained so far to a one- dimensional Bose fluid. Its density 
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operator is given by eq. (1) if we identify Qf /ft = po = pi 
(p = 2): 

f 

Psf = -d x ip + po (1 + cos(2(<^ + Qx))) + ... (69) 

ft 

d x ip is conjugate to the phase 9 of the Bose field op- 
erator—. Keeping our definitions of K, t and u, v de- 
notes now the phase velocity of the sound waves with 
v = \J n/(pom) and the elastic constant is c = n/(ftpo) 2 , 
where n is the compressibility per unit length (see also 
tabled). 



With the replacements 

K -> K~ l 
t -> t/K 2 
P = 2, 

(|18f) describes the action of the lD-supcrfluid in a random 
potential. The correlation functions for the superfluid 
can be obtained correspondingly from this replacements. 
To avoid confusions we write down the full action in this 
case explicitly: 



L K/t 

-F = ^E l dx I dT 



a,f3 



(70) 



K 



K/t 



S a B - — / dr' 



u cos 2 



(p a (x,T) - ipp(x,T')^j + <rd x ip a (x,T)d x (pp(x,T') 



Hence the RG-equations follow from (|19afl to l|19d|) 
with the above given replacements: 



dt 
~dl 
dK 

If 

du?_ 

~dT 

da 



1 4 

8u 2 



!6« „ ( llT , K\ , K 
— ttB 4/K, — coth — 
K 2 V ' 2t I 2t 



, K\ K 



2t 



2 K 

3 coth — 

K 2t 



supe 

calized phase occurs at K* = 2/3 Thermal fluctu- 
ations again suppress the disorder and destroy the the 
superfluid localization transition in ID. 



The model has q degenerate classical ground states 
given by tp m = 2irm/q with m = 0, q— 1. Performing a 
calculation analogous to the one above (but with u = 0) 
the RG-flow equations read 



dK 
~dT 
dt 
~dl 
dw 
~dJ 



„4„„2 D / „2 is \ CQth . 



2 ' V M rt 2t 



2r 



7T /i ■> „ / o if \ , K 
„4„„2 D . /„2^ \ CQth 



2f 



9 -if coth K 



4 2t 
where B\^ are given in H2U|I with 



to, 



9i 

!)2 



2x 2 cosx , 

(x 2 + t 2 ) cosx . 



(72a) 
t, (72b) 

(72c) 



VI. LATTICE POTENTIAL 

If the wave length A of the CDW modulation is com- 
mensurate with the period a (= ft, due to dimensionlcss 
units) of the underlying lattice such that nX = qa with 
integers n and q, the Umklapp term —2ft(w/K) cos(qip) 
appears in the Hamiltoniani. Therefore we switch on 
the lattice potential w ^ now. In this section we con- 
sider the case u = which leads to the Sine-Gordon type 
model: 



Slp 



L K/t 

dx I dr 



o 



2^ { i Q x f) 2 + (9^) 2 }-^ cos(^) 

(71) 



Plots of the functions B\ and B 2 can be found at the end 
of Appendix lAl 

At zero temperature <|72a|) and (|72c|) reduce to 



dK 

~dT 

dw 

~dT 



= -q i w 2 B 2 [q z K, oo), 



2 - —K 
4 



(73a) 
(73b) 



and we find, that for u = the lattice potential becomes 
relevant (i.e., w grows) for K < K w , where K w is implic- 
itly defined by 



w 2 (K u 



2ft q 2 f] \K* 



(74) 



12 



which follows from 



dK q 4 irijw \ J 

where we used l|73a|) and (|73b|) and the initial condition 
w(K* w = 8/q 2 ) = wo = 0; fj = -B 2 (q 2 K*, oo) (« 0.4, 
for g = 1). 

In this region the periodic potential stabilizes true 
long-range order of the CDW: the phase is everywhere 
close to one of the q classical ground states ip m . The 
depinning transition from the lattice for K — > K w — is 
again of KT type. The correlation length in the low- if 
ordered phase £ w is defined by w(]n£ w ) f=a 1 and diverges 
at K w — 0] . This can be seen by considerations anal- 
ogous to the disordered case. Defining 

7 = 2*1-2, 
K* 

2 77 6 ~ 2 

X = g9 rrw 

(note that fj > 0) leads for | | <C 1, i.e., close to K^, to 
the KT equations (|28b|) and (|25af) . 




FIG. 4: Typical flow diagram for the disorder free model in 
the three dimensional parameter space of K, w and t. w 
denotes the strength of the commensurate lattice potential 



At finite temperatures we find a similar scenario as 
in the case where we considered the influence of the 
disorder: w first increases in a K—t region given by 
M w = {(K,t)\K* > K cothf > 0}, i.e., when the r.h.s. 
of H72c|) is positive, but then decreases and flows into the 
region of large t and small w. Thus the periodic poten- 
tial becomes irrelevant at finite temperatures. This can 
be understood as follows: at finite t the ID quantum 
sine-Gordon model can be mapped on the Coulomb gas 
model on a torus of perimeter At since periodic bound- 
ary conditions apply now in the r-direction. Whereas 
the entropy of two opposite charges increases for sepa- 
ration L > At as ln(LAT), their action increases lin- 
early with L. Thus, the charges remain bound. The one- 
dimensional Coulomb gas has indeed only an insulating 
phase — . 



VII. CONCLUSION 

To conclude we have shown that in one-dimensional 
charge and spin density waves, Luttingcr liquids and 
suprafluids, quantum phase transitions between a dis- 
ordered (or locked-in) phase and an asymptotically free 
phase at zero temperature are destroyed by thermal fluc- 
tuations leaving behind a rich crossover behavior. This 
was demonstrated by using a full finite temperature 
rcnormalization group (RG) calculation. The crossover 
regions were characterized by the behavior of the phase 
pair correlation functions. For vanishing quantum fluctu- 
ations our calculation was amended by an exact solution 
in the case of strong disorder and by a mapping onto the 
Burgers equation with noise in the case of weak disorder, 
respectively. Both methods gave an exponential decay of 
density correlations. 

We have also briefly discussed, that the Coulomb in- 
teraction may destroy the unpinning (localization) tran- 
sition. 

The finite temperature calculation used in the present 
paper is also suited for treating the low frequency low 
temperature behavior of dynamical properties which may 
depend crucially on the ratio co/T. This as well as the 
discussion of the influence of quantum phase slipsi^ will 
be postponed to a forthcoming publication. 

The combined effect of disorder and the lattice poten- 
tial on the zero temperature phase diagram is still con- 
troversely discussed^^ and cannot be explained by the 
RG-results presented in this paper, since both pertuba- 
tions become relevant for small K. However, this prob- 
lem is beyond the scope of the present work and may be 
discussed in a future publication. 
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APPENDIX A: RENORMALIZATION OF THE 
DISORDER 

We present a short overview of the used finite temper- 
ature anisotropic renormalization group method in the 
case of the replicated disorder term. Starting point is 
the action 1)18(1. 

The phase field ip(x,r) = 

lrf3E^£| fc |<ie l( "" r+fer W„ with the Matsub- 
ara frequencies ui n = 2imt/K and k = 27rm/(LA) (note 
that rescaled coordinates are used) is split in a slow 
(|fc| < b" 1 ) and a fast mode part (6" 1 < \k\ < 1), where 
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b = e is a rcscaling parameter of order 1 . Notice, that 
the y?< still have all Matsubara Fourier components. 



In order to find the RG-corrections of the other pa- 
rameters in the model, we follow Wilson^ and expand 




where R[f] = cos(p/). ((. . .}) > denotes the cumulative 



or connected average over the fast modes in the " momen- 
tum stripes" with the free gaussian model. The correc- 
tion in first order is given by 




For calculating the cumulative average of the func- 
tional R, R is expanded in small Aip > = tp ai> (x,T) — 
Pft>(iy). e.g. 



{(R[A<p])) 0>> = -p 2 {R[A v< ] <^ >> > 0)> -RV*p<] ^a,>(x,T)^,>( aJ y)) 0)> ) +0(A<p*). (A3) 

The first term in (S^J gives a correction to the disorder parameter u and the second term a correction to K. The 
free thermal average over the fast modes can be evaluated with the free propagator (k 2 +u>f l )~ 1 and using the formula 



cos(nx) 7r cosh((-7r — x)a) 
n 2 + a 2 | a | : 

to treat the sum over the Matsubara frequencies yielding 



Ecos(nx) 7r cosh 7r-x a 
n 2 + a 2 \a\ smh(7r|a|) 



n— — oo 



, , m K cosh(K/2t - \Ar\j . , , 



with At = t — t'. 

In order to find a good gaussian approximation for R, we perform a variational calculation for a Sine Gordon model: 

Use = n a + Hi = Ho - n J d d r cos(p(r)) , (A5) 
where Ho is the gaussian part. We approximate H\ by 

Hl = |^ r ^ 2(r)j (A6) 

and define H = Ho + Hi- To find the optimal function «(r), the variational free energy}^ J^var, with 



Fsg < T var = F+(H SG -H) (A7) 



is minimized with respect to k: 



= | Jd d r («(r) - /icosMr)))^ «¥>W(f)»* . (A8) 

For the last equality we took into account that the averages are gaussian such that we could apply the Wick theorem. 
From (|A8|> we finally get 

«(r) = n (cos(¥>(r)))^ = ^cxp (~ (ip 2 (r))^j . (A9) 
For small disorder (|A9|) yields for i?[Ayj<] the approximate expression 

i?[A^<] ~ -^(A^) 2 e -^< (Av < )2 )o- . (A10) 
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The same result can be obtained in terms of an operator product expansion^ of R[f]. In order to get a RG correction 
to K a gradient expansion of (Ay> < ) 2 in l|A10|) in small At is performed, which is justified by the exponential decay 
of the correlation function (|A4() on the integration interval such that higher orders in At do not contribute to the RG 

correction: (A^<) 2 « (d T <p<(x, T)At) 2 with T = (t + t')/2. 

The pair correlation function in the argument of the exponential function in IjAlOf) can be approximated by the 
expression already shown in (|55|l . i.e., 



<0( 



X,T 



^(o,o)) 2 ) f 



K 



In 1 



(-) 

\2irt J 



2 r 



cosh 



[IT) 



(All) 



After integration of the fast modes of tp, one rescales 
the system to maintain the fluctuation strength and the 
spatial density of the degrees of freedom (I = \nb): 

x — > x' = xb~ x 

t — > t' = rb- z , T ^ T' =Tb z 
tp — > tp' — ipb" 1 ' 

which leads to rcscalcd parameters. For our model these 
are given by (here in d dimensions) 



c 


= cb d+z - 2+2 ^ 


(A12a) 


V 


= vb z ~\ 


(A12b) 


A" 


= Kb 1 -^, 


(A12c) 


t' 


= tb 2 - d - 2 ^ 


(A12d) 


u 


= ub 2 - d / 2 , 


(Al2e) 


a' 


= ab 2 - d , 


(Al2f) 


w 


= wb 2d . 


(Al2g) 



Due to the invariance of the system under a phase shift 
of 2mr, £ is zero for symmetry reasons. 

The RG contribution to the flow equation for u 2 follows 
from the first term of (|A3|) and l|A4|) with At = 0: 



u 2 [ 1 



coth — In b 

2 2t 



Together with (|A12e|) (d = 1) one gets 



dl 



p 2 K , A 

coth — 

2 2t 



(A13) 



(A14) 



The RG correction to K follows from the second term 

of with inn and mni: 



A' = A 1- 



P 



K/(2t) 

a , K f , 9 cosh(f - 

-u 2 coth— / dTT 2 

2 2t J cosh K 





2/ 



e -^((^(0,r)- V (0,0)) 2 ) o l„fcj f (A15) 

and the flow equation (for d = 1 + e) follows from l|A12cf) : 



A' - A 



A 



A 



-e - t" 2 coth ^-Bo ( p 2 K 

2 2t V 2i 



A 

(A16) 




FIG. 5: Function Bo(p 2 K, ||) plotted with respect to K for 
different temperatures. 



with Bq given in (|20|) . for which wc used IjAlip . 

The function Bq which appears in this flow equation is 
plotted in Fig. [S] 

For completeness we also plot the functions Bi and B 2 
in the relevant A-region for the lattice unpinning transi- 
tion [see eqs. (I72all and (|72b|) ]. Note, that for evaluating 
these functions at zero temperature, one has to execute 
the integrals at finite temperature first and then take the 
limit t -> 0. 



APPENDIX B: STRONG PINNING 

To calculate the phase correlation function in the 
strong pinning limit, it is necessary to study the order 
statistics of the impurity distances ti — {xi—Xi-\) — li mp . 
Following David2£, we obtain for the probability den- 
sity function (pdf) of the e^'s in the case of uniformly 
distributed impurity positions < xt < L : p(ei) = 

Kmpi 1 - T( e i + h,np)) N ~ 1 with -l imp < Q < 

the thermodynamic limit the pdf can be rewritten as 



' imp ■ In 



p(ei) ~ -^-e l """ ei , -l imp < £i < oo . (Bl) 
e 

With this, one can calculate the n-th moment e™ (for 
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B-|(q 2 K,K/2t), q=1 
2 4 6 8 10 12 14 




B 2 (q2K,K/2t), q=1 
2 4 6 




FIG. 6: Function B 1 (q 2 K, §) and B 2 (q 2 K, §) plotted with 
respect to K for different temperatures, written next to the 
graphs, and q — 1. 



n > 1, ej = 0) as follows: 



kmp 6 

e 1 J <9A™ 



(B2) 



yields: e™ = ^ X)fc=2 ^TT - an d f° r tnc correlator eie^ 



I 2 6- 



With this results we can derive the pair correlation 
function i|tj3[l • Therefore we calculate the discrete version 



(ifn — tpi) ) in the limit T —> 0. With the definitions 



given below eq. I|3b|) . we can rewrite 

2 /n-l 



2 4t^ ('^ f, pQe t 



27T 



(B3) 



Using e; = /ij = and hihj oc <5y leads to 



((^+i-^) 2 ) = ^(V^ 



7 (B4) 



=c 



Because only the value 



2tt t ' 



for hi is taken into 



account for evaluation of the thermal average, we get 



dc, 



2tt 



l imp ( pQa , 

e rmp h 7 — 

/. e V 2-7T 



2ir 



+ 7 



If we substitute x = l ir ^ p ei + 1 and take into account that 
[x + n] = [x] + n for n G Z, we get 



C = / da; e 
./o 



V2a 



2a 



(B5) 



with the parameter a = — . Now the quadratic term 
in the integral is expanded, which leads to the following 
three (converging) integrals 



h 
h 



i r. _ l2 i 

— T / dx e x x z = — T 
4a 2 J 2a 2 



dx e x x 

(2fc+l)a 



X 

2T\ 



E 

fc=i 



dx e x x 



x 

2a~ 



(B6) 



(2fc-l)a 



E 

k=l 



dxe x 

(2fe+l)< 

/ 

(2fc-l)a 



2a 





' X ' 


y dxe x 






2a- 



with C = 7i +/ 2 + I 3 . For (2fc - l)a < a: < (2k + l)a, 
fceZ. [^-1 = fc, such that the Gaussian brackets in the 
finite integrals in /2 and ^3 can be replaced by k or k 2 , 
respectively: 



I 2 = 

*=i 
(2A 

00 „ 

/ 3 = E^ 2 / 

7 1 «* 



(2fc+l)a 

dx e~ x x , 

(2fc-l)a 
(2fc+l)a 

dx e~ x . 



(B7) 



k=l 



(2fc-l)a 

The values of these two simple integrals are 

(2fc+l)a 



J dxe~ x x = 2e~ 2fea ((l + 2/sa)sinh(a) 
a cosh(a)) , 
= 2e~ 2fca sinh(a) . 



(2fc-l)a 
(2fe+l)a 



dxe 



(2fc-l)a 

The remaining sums in I2 and ^3 are only derivatives 
of the geometric scries which can be easily evaluated. 
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The result is 



C 



1 \ 



2a \ a sinh(a) 



(B8) 



1 

12 



7 

720' 



31 



30240 



0{a 6 



where the expansion is useful only if a is small, i.e., 
QUmp 3> 1- If we would have neglected the order statis- 
tics of the impurity distances we would get only the lead- 
ing constant: C = 1/12. 

I|B8(I yields the presented expression for the pair cor- 
relation function JHSJl- 



APPENDIX C: CORRELATION LENGTH IN 
THE CLASSICAL REGION AT FINITE 
TEMPERATURE 

In the weak pinning limit, 1 <C hmpQ <C c/(Ui mp pi) or 
Lfl ^ Ump, the classical Hamiltonian can be rewritten 
to a random-field XY-model 



HdassiL) = / dx < | {d x <p(x) - a) 2 + 



U imp Pi 



cos\p(<p(x) - a(x))} (CI) 



imp 



where a(x) is a random phase with zero average and 



sV (a(*)-a(»0) = t g( x _ x >y 



(C2) 



In the following we consider only the backward scattering 
term of the correlation function Cj, and therefore neglect 
the forward scattering term a = Ui mp f / (cnli mp ). 

The goal is now to find an effective temperature, which 
replaces T in the correlation function for the free case 
(I64|l . We start with a Burgers-like equation, which one 
gets after a Cole-Hopf transformation from the transfer 
matrix equation, for the restricted free energy J-(x, ip) = 



-T hi Z(x,(p) with the partition function 

tp(x)=ip 



Z(x,<p)= J Vtpt 



-H(x)/T 



V(0) 



The equation reads 



dT T d 2 T 1 (dTY U, 



dx 2c d(p 2 2c \ dip J I 



-— ( — ) H — cos[p(tp{x) - a{x))} 



l nnp 



U(.x,tp) 



Using the Fourier transform F(x, ip) 



(C3) 



dk duj i(uj<p — kx) 



f dk c 
J (2tt 



rewritten as 



T{k,u) (analogous for U(x,pj), l|C3|l is 



-ikT{k,u) = -^T{k,ui) + U(k,u)+ (C4) 
/ <g$V(w - io')T(k - k', lo - u')F{k', to') , 



with 



U{k,u) 



h+{k)5(uj — p) 



limp \ 

+h_{k)5(uj + p) \, (C5) 



Introduction the dimensionless quantities 



irtu> 2 /2 - ik/K 

U imp Pi 



i(k,u>) = (eAc)~ 1 C/(/c,w), 



(C6) 

(C7) 
(C8) 



and setting J-{k,uj) = ceg(k 7 u>)u(k,uj) we obtain the fol- 
lowing, self-consistent equation for the Green's function 



g(k,uj)u(k,uj) = g {k,uj)u(k,Lu) + ^g Q (k,uj) J ^^~-^u'(uj-u')g(k-k',u-u')g(k',u}')u(k-k' ,u-u')u(k' , J) (C9) 

For e < 1, i.e., for weak disorder, this equation is iterated to first non-vanishing order in e {one-loop approximation) 
and averaged over disorder. The disorder average u(k,u>)u(k' ,oj') can be calculated using (|C5|I and l|C2|l . which gives 

u(k,u))u(k\uj>) = AV (h + (k)h-(k')S(u} - p)5(u' + p) + h-{k)h + {k')5{uj + p)5{u' - p)\ 
= 2A 2 l„ np 'K 3 5(k + k')S(uj + u)') {5(uj + p) + S(uj - p)} 

= 2S(k + k')S(w + u')D(w, k) . (CIO) 
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Therefore, we get for g in order e 2 



l(k,w) =g {k,u) + e g (k,u>) 



dk' dw 1 
A 2 (2tt) z 



{ui-u;'yu;(-u')go{k\w')go(k- k\w-u;')g {-k',-u;')D(k',u;'). (Cll) 



fa; 




g(k,a) 



j U(/C-/c',(B-(0') 

g(k-k',ay-(s!) 

g(k',m') 
c-j u(k',(o') 

D{kW) 



+ 2 



/Ox 



/(,C0 /C-/C',(B-(B' /c,co 

FIG. 7: Diagrams for (a) eq. £U and (b) eq. iCTfl . 

The diagrams visualizing eqs. (|C9j) and I|C11|1 are depicted in figured 
For k = (|C11|I reduces to 



KO^) = 5 o(0,w) 



1 + e 



2 P%np J 

87rfcj 



1 



p — u> 



p + uj 



(fc'/A) 2 + (ip 2 i/2) 2 [7rt(w-p) 2 /2-ifc'/A 7rf(w + p) 2 /2 - zfc'/A 
I 



(C12) 



Because we calculate the correlation length in the ther- 
mal regime (see figure [3| with t > t u the fc'-integral in 
(|CJ1 21 gives the biggest contribution to g at small u>. In 
this case the /c'-integral can be easily calculated which 
leads to 



fl(0,w«l)« fl o(O l «) 1 



2 A.1>imp 



2i: 3 p 2 t 3 



(C13) 



or for the effective Temperature T eJ y 
1 



1 f\. ^2 ^-^imp 



2 \7rp 2 t 



which yields for the correlation length 



(C14) 



c 1 « 2 /(T) * 



2 VttP 2 * 



as written in the text. For high temperatures t S> t u we 
recover the linear t dependency of the free case. 

A related calculation for directed polymers and inter- 
face growth can be found in j3^. Note, that in this paper 
x plays the role of <p and t the role of x in the above cal- 
culation. 



APPENDIX D: SYMBOL REFERENCE 



A, (C15) 
TABLE III: List of used quantities 



symbol 


quantity 


eq. ref. 


a 


lattice constant 




b = e~' 


rescaling parameter of order 1 




Bi 


functions used in flow equations 


CP 


c 


elastic constant 


teat 


C(x,r), C f , C b 


phase correlation functions 


0, C23 
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TABLE III: (continued...) 



symbol 


quantity 


eq. ref. 


d 


spatial dimension 




f(T) 


condensate density 


i—i 


■rr 


free energy 




9 


dimensionless electron-phonon coupling constant 


rn 

w 


9i 


functions used in Bi 




<hi 


electron-phonon coupling constants 




hi 


integers (cf. strong pinning) 


(O 


~H, Ho 


Hamiltonian (complete and free) 


m 




wave vectors 




k.F 


Fermi wave vector 




K, K(l), K 


dimensionless parameter for quantum fluctuations 


(fBbfl . 1101)0 


K u , K* 


K- values denning the separatrix/fixed point of the disorder unpinning transition 






K- values defining the separatrix/fixed point of the lattice unpinning transition 


cm 


L 


system length 


Q 


Lfl 


Fukuyama-Lee length 





limp 


mean impurity distance 




r 

Lo 


free (gaussian) part of the Lagrangian 




Hi 


integers (cf. strong pinning) 


till 


AT 

1 * imp 


number of impurities 


■1 1 '1 


P 


commensurability used in the density 


Ul 


n(fA 


T^rnriPi ni 1 \^~\t HpnQit^r fuTirfinti cit 




r 


momentum operator, conjugate to (p 


,1 ( | 


Q 


density wave vector 


W 


<7 


commensurability used in the lattice potential 


(I4cl) 


o o 
<->, oo 


action (full and gaussian part) 




S y ' 


replicated action 


d 


Ssf, ShP 


action for superfluids and lattice potential, respectively 


czni 


S, Si 


density correlation functions 




T 


temperature 






mean-field condensation temperature 




t, t(l), t 


parameter for thermal fluctuations 


Eat 


t u = 1/(ALfl) 


crossover temperature from classical disordered to thermal regime 




t K 


temperature separating the thermal and disordered regime 




U(x) 


disorder potential 


ph)l 


Ui 


impurity potential 




Uimp 


mean impurity potential 




u, u(l), u 


dimensionless parameter for disorder fluctuations 




I'F 


Fermi velocity 




V 


phason velocity 





V c (x) 


Coulomb potential 


ca 


W 


lattice potential strength 


1 r 


w, w(l), Wo 


dimensionless parameter for lattice potential strength 


JSHl 




impurity positions 




z 


dimensionless distance in r-x-space 




(X 

(3 


parameter used in the strong pinning limit 
inverse temperature 




7 


parameter for KT flow equations 


if^al 


A 


order parameter 


in 


ti 


deviation from mean impurity distance 


EHl 


C 


transverse width of the quasi one-dimensional system 




V, V 


= Bo(p 2 K*,oo), = —B2(q 2 /K^,co), respectively 




A 


density wave length 




At 


de Broglie wave length 
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TABLE III: (continued...) 

symbol quantity eq. ref. 

A momentum cutoff 

£, £ U) £ w , etc. correlation lengths tabic Iffl 

p(x), psf{x) charge/spin or superfluid density Q, 16911 

po mean density 

Pi density amplitude for harmonic part of p(x) Q 

a forward scattering amplitude 1181 

r imaginary time coordinate 

ip phase variable @ 

X parameter for KT flow equations l|28bfl 

T auxiliary function 12111 

uj n Matsubara frequencies 
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